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PARALLELIZATION OF MODULAR ALGORITHMS 



NAZERAN IDREES, GERHARD PFISTER, AND STEFAN STEIDEL 



Abstract. In this paper we investigate the parallelization of two modular 
algorithms. In fact, we consider the modular computation of Grobner bases 
(resp. standard bases) and the modular computation of the associated primes 
of a zero— dimensional ideal and describe their parallel implementation in Sin- 
gular. Our modular algorithms to solve problems over Q mainly consist of 
three parts, solving the problem modulo p for several primes p, lifting the result 
to Q by applying Chinese remainder resp. rational reconstruction, and a part 
of verification. Arnold proved using the Hilbert function that the verification 

■^^ part in the modular algorithm to compute Grobner bases can be simplified 

for homogeneous ideals (cf. A03 ). The idea of the proof could easily be 

^_ ■ adapted to the local case, i.e. for local orderings and not necessarily homo- 

geneous ideals, using the Hilbert-Samuel function (cf. |Pf07l ). In this paper 
we prove the corresponding theorem for non— homogeneous ideals in case of a 
global ordering. 



k*" ' 1. Introduction 

m : 

^»0 , We consider an ideal in a polynomial ring over the rationals. In section [2] we 

\& ' describe a parallel modular implementation of the Grobner basis (resp. standard 

basis) algorithm. Afterwards we restrict ourselves to the case of a zero-dimensional 
l/"~) ' ideal and introduce a parallel modular implementation of the algorithm to compute 

the associated primes in section [3] Finally we give a couple of examples with 
corresponding timings and some conclusions in section |4] Both algorithms are 
implemented in Singular. The Grobner basis resp. standard basis algorithm can 
be found in the library modstd . lib and the algorithm for computing the associated 
primes in assprimeszerodim.lib. They are included in the release Singular 3- 



U ' 1-2. 



The task to compute a Grobner basis G of an ideal I using modular methods 
consists of three steps. In the first step, we compute the Grobner basis modulo p 
for sufficiently many primes p and, in the second step, use Chinese remainder and 
rational reconstruction to obtain a result over Q. In the third step, we have to 
verify that the result obtained this way is correct, i.e. to verify that / = (G) and G 
is a Grobner basis of (G). If this fails we go back to the first step. The third step is 
usually at least as time consuming as the first step. Omitting the third step would 
produce a Grobner basis only with high probability and the result could be wrong 
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in extreme situations. It is known that some of the commercial computer algebra 
systems have problems in this directions 

Arnold proved using the Hilbert function that the verification part in the mod- 
ular algorithm to compute Grobner bases can be simplified for homogeneous ideals 
(cf. [A03 ): Let / C Q[xi, . . . , x n ] be a homogeneous ideal, > a global monomial 
ordering and G C Q[x\, . . . , x n ] be a set of polynomials such that / C (G), G is a 
Grobner basis of (G) and LM(G) = LM(I¥ p [xi, . . . , x n ]) for some prime number 
p where LM(G) denotes the set of leading monomials of G w.r.t. >, then G is 
a Grobner basis of /. The idea of the proof could easily be adapted to the local 
case, i.e. for local orderings and / not necessarily homogeneous, using the Hilbert- 
Samuel function (cf. [Pf07] ). In this paper we prove the corresponding theorem 
for non-homogeneous ideals in case of a global ordering. Two important assump- 
tions of the theorem are the facts that I C. (G) and G is a Grobner basis of (G). 
This verification can be very time consuming in a negative case. Hence, we use 
a so-called pTestSB which is one of the new ideas for our algorithm. Therefore 
we randomly choose a prime number p which has not been used in the previous 
computations and perform the verification modulo p. Only if the pTestSB is 
positive we perform the verification over Q, and the last required condition that 
LM(G) = LM(I¥ p [xi, . . . , x n ]) is then automatically fulfilled. 

The implementation of our algorithm as Singular library implies that we did 
not change the kernel routines of Singular. We plan to implement the algorithm 
in the kernel of Singular in future. For this purpose we can apply the ideas of 
Grabe (cf. |G94| ) - using multimodular coefficients - and Traverso (cf. JT89] ) - 
using the trace-algorithm. The trace-algorithm would speed up the computations 
in positive characteristic a lot. We compute a Grobner basis of an ideal I C 
Q[xi, . . • , x n ] over ¥ p [xi, . . . , x n ] for a random prime p and keep in mind the zero- 
reductions of the s-polynomials such that we do not perform these reductions in 
any other Grobner basis computation over ¥ q [xi, . . . , x n ] for primes q ^ p. We do 
not need these information, i.e. the guarantee that we really obtain a Grobner basis 
over ¥ q [xi, . . . , x n ], since we have the verification step - that the lifted result over 
Q[xi, . . . , x n ] is a Grobner basis of I - at the end anyway. 

Our idea regarding the primary decomposition of a zero-dimensional ideal I C 
Q[X] is to compute the associated primes M\,...,M S of I and use seperators 
<7i, . . . , crjj such that the saturation / : <jf° of / w.r.t. Oi is the primary ideal corre- 
sponding to Mi (cf. [SY96] ). The computation of the associated primes is based on 
the so-called Shape Lemma (Proposition[3jj2)). Here, one new idea is to choose a 
generic linear form r — a\X\ + . . . + a n —iX n —i +x n with a\, . . . , a n — l € Z and a ran- 
dom prime p to test if dim Fp (F p [X] /I¥ p [X]) = dim Fp (¥ p [x n ]/(ip(I)¥ p [X]Pi¥ p [Xn])), 
i.e. i/j(I)¥ p [X] = (x\ — gx(x n ), . . . , X n -% — g n -i(x n ), F(x n )) whereat ip denotes the 
linear map defined by ip{xi) = Xi for i = 1, . . . , n—1 and ip(x n ) = 2x n — r. If this test 



Let N be the product of all primes smaller than 2 32 and I = (v + w + x + y + z, vw + wx + 
xy + yz + vz, vwx + wxy + xyz + vyz + vwz, vwxy + wxyz + vxyz + vwyz + vwxz, vwxyz + TV) C 
Q[v,w,x,y,z}. Then MAGMA V2. 16-11 (64-bit version) computes a wrong Grobner basis, in 
particular it computes the Grobner basis of the ideal J = (v + w + x + y + z, vw + wx + xy + yz + 
vz, vwx + wxy + xyz + vyz + vwz, vwxy + wxyz + vxyz + vwyz + vwxz, vwxyz) C Q[v, w, x, y, z] 
which obviously differs from /. 

We call ai a seperator w.r.t. M; if a^ £ Mi and (X; £ Mj for j ^ i. 
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called pTestRad is positive then the ideal / in Q[X] has the same property with 
high probability. If the test is negative then we compute the radical of / using the 
idea of Krick and Logar fProposition l3.3F l)) combined with modular methods, and 
replace / by \fl. Afterwards we compute (F) = (I,T— r)q[x,T] H Q[T], again using 
modular methods, i.e. we compute F (p *> such that (F&) = (I,T-r) ¥p[x ,T] nF p [T] 
and deg(F^ p ') — dimq(Q[X]/I) for sufficiently many primes p, and we use Chi- 
nese remainder and rational reconstruction to obtain F £ Q[T]. The verification 
is the test whether it holds F(r) £ I and no proper factor of F(r) is in /. If 
F = Fp ■ ■ ■ F" s is the factorization of F in Q[T] into irreducible factors then 
M i = (I, F\ (r)} , . . . , M s = (I,F S (r)) are the associated primes of /. The new ideas 
in this approach are the pTestRad described above and the fact that we do not 
compute the associated primes in positive characteristic but instead one special 
generator of the radical, F(r), which is much better to controljj 

We use the following notation. Let X = {xi, . . . , x n } be a set of variables. We 
denote by Mon(X) the set of monomials, and by Q[X] the polynomial ring over Q 
in these n indeterminates. Let S C Q[X] be a set of polynomials, then LM(5) := 
{LM(/) | / € S} is the set of leading monomials of S. Given an ideal / C Q[X] 
we can always choose a finite set of polynomials Fj such that / = (Fj). li I = 
(/i , . . . , f r ) C Q[X] and p is a prime number which does not divide any denominator 
of the coefficients of fi, ■ ■ ■ , f r we will write I p :— (/i mod p, . . . , f r mod p) C 
V P [X]. 

2. Computing Grobner bases using modular methods 

In the following we consider an ideal I = (/i, . . • , f r ) C Q[X] together with a 
monomial ordering > and set Fj = {/i, . . . , f r }. We assume that > is either global 
or local. Within this section we describe an algorithm for computing a Grobner 
basis resp. a standard basif|j G C Q[X] of / by using modular methods. 

The basic idea of the algorithm is as follows. Choose a set P of prime numbers, 
compute standard bases G p of I p C F P LY], for every p £ P, and finally lift these 
modular standard bases to a standard basis G C Q[X] of /. The lifting process 
consists of two steps. Firstly, the set GP := {G p | p £ P} is lifted to Gn Q 
Z/NZ[X] with N := Yi P epP by applying the Chinese remainder algorithm to the 
coefficients of the polynomials occuring in GP. Since Gn is uniquely determined 
modulo TV, theory requires N to be larger than the moduli of all coefficients occuring 
in a standard basis of I over Q. This issue is not revisable a priori and will be 
discussed later in this section. Secondly, we obtain G C Q[X] by pulling back the 
modular coefficients occuring in Gn to rational coefficients via the Farey rational 
mapj. This map is guaranteed to be bijective provided that \/N/2 is larger than 
the moduli of all coefficients in GO The latter condition on N concerning the 
Farey rational map obviously implies the former condition concerning the Chinese 



The computation of the associated primes in positive characteristic would create similar prob- 
lems as the factorization of polynomials: Different behaviour of splitting in different characteristics. 
Therefore it is easier and faster to compute F 6 Q[T] and factorize this polynomial. 

For definitions and properties cf. [GP07 . 
5 Farey fractions refer to rational reconstruction. A definition of Farey fractions and the Farey 
rational map can be found in A03 , KG83 , Pf07); for remarks concerning its computation cf. 
|KG83| . 

Remarks on the required bound on the coefficients are given in |KG83| . 
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remainder algorithm. We consequently define two corresponding notions that are 
essential regarding the algorithm. 

Definition 2.1. Let G be a standard basis of I. 

(1) f f G p is a standard basis of I p , then the prime number p is called lucky for 
I if and only if LM(G) = LM(G P ). Otherwise p is called unlucky for I. 

(2) A set P of lucky primes for I is called sufficiently large for I if and only if 
H„ e pP > max{2 • |c| 2 | c coefficient occuring in G}. 

Now we can concretize the theoretical idea of the algorithm. Consider a suffi- 
ciently large set P of lucky primes for / such that none of these primes divides any 
coefficient occuring in Fj, compute the set GP, and lift this result to a rational 
standard basis G of / as aforementioned. More details can be found in [A03 . 

In practice, we have to handle two difficulties since naturally the standard basis 
G of I is a priori unknown. In fact, it is necessary to ensure that every prime number 
used is lucky for /, and to decide whether the chosen set of primes is sufficiently 
large for /. 

Therefore, we fix a natural number s and an arbitrary set of primes P of cardi- 
nality s. After having computed the set of standard bases GP := {G p | p £ P} we 
delete the unlucky primes in the following way. 

deleteUnluckyPrimesSB: We define an equivalence relation on (GP,P) by 
(G p ,p) ~ (G q ,q) :<=> LM(Gp) = LM(G ? ). Then the equivalence class of largest 
cardinality is stored in (GP,P), the others are deleted. 

With the aid of this method we are able to choose a set of lucky primes with 
high probability. A faulty decision will be compensated by subsequent tests. 

Since we cannot predict if a given set of primes P is sufficiently large for /, we 
have to proceed by trial and error. Hence, we lift the set GP to G C Q[A], as 
per the description at the beginning of this section, and test whether G is already 
a standard basis of I. Otherwise we enlarge the set P by s new prime numbers 
and continue analogously until once the test is positive. The test especially verifies 
whether G is a standard basis of (G), but this computation in Q[X] can be very 
expensive if P is far away from being sufficiently large for /. Hence, we prefix a test 
in positive characteristic that is a sufficient criterion if P is not sufficiently large 
for I. 

pTestSB: We randomly choose a prime number p £ P such that p does not 
divide the numerator and denominator of any coefficient occuring in Fj . The test 
is positive if and only if (G mod p) is a standard basis of I p . We explicitly test 
whether (fi mod p) £ (G mod p) for i — 1, . . . ,r and (G mod p) C std(/ p jj. 

This test in positive characteristic accelerates the algorithm enormously. It is 
much faster than in characteristic zero since the standard basis computation in 
pTestSB is as expensive as in any other positive characteristic, i.e., as any other 
standard basis computation within the algorithm. 

If the pTestSB is negative, then P is not sufficiently large for /, that is, G 
cannot be a standard basis of / over Q. Contrariwise, if the pTestSB is positive, 
then G is most probably a standard basis of /. 

The procedure std is implemented in Singular and computes a Grobner basis resp. standard 
basis of the input. 
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Algorithm Q] shows the modular standard basis algorithm^ 

Algorithm 1 modStd 

Assume that > is either a global or a local monomial ordering. 
Input: ICQ[X]. 
Output: G C Q[X] the standard basis of /. 

choose P, a list of random primes; 

GP = %- 

loop 

for pePdo 

compute a standard basis G p of I p ; 
GP = GPU{G P }; 
(GP, P) = deleteUnluckyPrimesSB(GP, P); 

lift (GP,P) to G C Q[A] by applying Chinese remainder and Farey rational 
map; 

if pTestSB(7, G,P) then 
if I C (G) then 

if G is a standard basis of (G) then 
return G; 
enlarge P; 



Remark 2.2. The presented version of the algorithm is just pseudo-code whereas 
its implementation in Singular is optimized. E.g., the standard bases G p of 
I p C F P [X] for p 6 P arc not computed repeatedly, but stored and reused in 
further iteration steps. 

Remark 2.3. Algorithm Q] can easily be parallelized in the following way: 

(1) Compute the standard bases G p in parallel. 

(2) Parallelize the final tests: 

• Check if I C (G) by checking if / e (G) for all / € Fi. 

• Check if G is a standard basis of (G) by checking if every s-polynomial 
not excluded by well-known criteria, vanishes by reduction w.r.t. G. 

Algorithm [T] terminates by construction, and its correctness is guaranteed by the 
following theorem which is proven in [A03] in the case that / is homogeneous resp. 
in [Pf07] in the case that the ordering is local. The case that the ordering is global 
follows by using weighted homogenization as in Theorem 7.5.1 of [GP07J . 

Theorem 2.4. Let G C Q[X] be a set of polynomials such that LM(G) = LM(G P ) 
where G v is a standard basis of I p for some prime number p, G is a standard basis 
of (G) and IC(G). Then I = (G) . 

Note that the first condition follows from a positive result of pTestSB whereas 
the second and third condition are verified explicitly at the end of the algorithm. 



Proof of Theorem \2.4\ We assume that > is a global monomial ordering. The proof 
for a local ordering is similar. Let Fj — {/i, . . . , f r } C Q[X] such that / = (Fi) 



°The corresponding procedures are implemented in Singular in the library modstd.lib. 
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and G = {gi, . . . ,g s } C Q[X]. Since G is a standard basis of (G) w.r.t. > and 
I C (G) there exist for each i = 1, . . . ,r polynomials fy G Q[-X"] such that 

s 

fi=^2 ^'* satis fy in g LM > (/t) > LM > teift) for all j = 1, . . . , s. 

Due to Corollary 1.7.9 of [GP07] there exists a finite set M C Mon(X) with the 
following property: Let >' be any monomial ordering on Mon(X) coinciding with 
> on M, then LM>(G) = LM>/(G) and G is also a standard basis of (G) w.r.t. >'. 
Moreover, due to Lemma 1.2.11 resp. Exercise 1.7.17 of GP07] we possibly 
enlarge the set M and chosse some w = (wi, . . . , w n ) £ Z™ such that >=>„, on 
M, i.e. LM>(G) = LM >ra (G) resp. G is a standard basis of (G) w.r.t. > w , andj 

w-deg(LM>„(/,)) > w-deg(LM >m (tail(/ i ))) , 

w-deg(LM >ra ( 5 ,)) > w-deg(LM >iu (tai%))) , 

w-deg(LM>„ feiSi)) > w-deg(LM >TO (tail(foft))) , 

for all ? = 1, . . . ,r and j — 1, . . . , s. 

Now we consider on Q[X, t] the weighted degree ordering with weight vector 
(ioi, . . . , «i„, 1) refined by >„, on QLY] and denote it also by > w . For / G Q[X] 
let f h = i w - dc «(/) • f{xx/t Wx ,...,x n /t Wn ) be the weighted homogenization of / 
w.r.t. t. We set Fj := {/ x ft , . . . , ft }, 7 := <F/> and G := {g^ . . . , g' s 1 } . Then 
Proposition 7.5.3 of |GP07) guarantees that G is a standard basis of (G) and since 
LM >u: (G) = LM >ra (G p ) it also holds by construction that LM >ro (G) = LM >tu (G p ). 
Now let i £ {1, . . . , r}, then /j = £^ =1 ?ij5j satisfying LM>„(/j) > w LM>„ (£#&) 
for all j = 1, . . . , s. This implies w-deg(/j) > w-deg^ygj) for all j = 1, . . . , s by 
the choice of w £ Z> . Consequently we have 

f w-dcg(/i) t (EL Eh 



t" 



3=1 



f w-dc g (/ l)e (EL ,..,,EE) g (EL^.^EE] e (g) 



thus Z/ 1 G (G) resp. J C (G) since ie {1, . . . , r} was arbitrarily chosen. 

It remains to prove that / = (G) . Let n £ N. We know that i p = (G p ) due to the 
fact that LM >ro (G) = LM >ra (G p ), so especially it holds HF 7p (n) = HF/^ \{n) = 
HF/gV (n) for the corresponding Hilbert functions. On the other hand we have 

HF 7 (n) < HF 7p = HF (Z?) (n) < HF» < oo, 

where the second inequality is true since I C (G). The first inequality follows from 
the fact that dimQ(J[n]) > dimw p (I p [n]), where I[n] resp. I P [n] denotes the vector 
space generated by all (weighted) homogeneous polynomials of degree n. Namely 
we can find a Q-basis of I[n] of polynomials in Z[X, t]P\I which induces generators 
of J p [n]. D 

Remark 2.5. Algorithm [T] is also applicable without applying the final tests, i.e. 
skipping the verification that / C (G) and G is a standard basis of (G) . In this case 
the algorithm is probabilistic, i.e. the output G is a standard basis of the input 



9 For a polynomial / G Q[X], we define by tail(/) := / - LM(/) the tail of /; cf. IGP07I . 
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/ only with high probability. This usually accelerates the algorithm enormously. 
Note that the probabilistic algorithm works for any ordering, i.e. also for the 
so-called mixed ordering. In case of a mixed ordering one could homogenize the 
ideal /, compute a standard basis using modStd and dehomogenize afterwards. 
Experiments showed that this is usually not efficient since the standard basis of the 
homogenized input has often much more elements than the standard basis of the 
ideal we started with. 

3. A MODULAR APPROACH TO PRIMARY DECOMPOSITION 

In the following let / C Q[X] be a zero-dimensional ideal and d :— dimq(Q[X] / 1) . 
Within this section we describe an algorithm for computing the associated primes 
of / using modular methods. In conclusion we make remarks how to achieve the 
corresponding primary ideals from the associated primes of I. 

The following well-known proposition (cf. GTZ88 or GP07 ) describes how 
to compute the associated prime ideals of a radical ideal over Q. Note that these 
results are also valid for perfect infinite fields. 



Proposition 3.1. Let I C Q[X] be a radical ideal. 

(1) Let (F) = In Q[x„] and assume deg(F) = dim Q (Q[X] / 1) . Let F = 
Fi ■ ■ ■ F s be the factorization of F into irreducible factors over Q. Then 
I = n"=i (I) ^i) an d (!■> Pi) * s prime for i — 1, . . . , s. 

(2) There exists a non-empty Zariski open subset U C Q" _1 such that for all 
a = (ax, . . . , a„_i) G U the linear coordinate change ip a defined by (f a (xi) = 
Xi for i = 1, . . . ,n — 1 and (f a {x n ) = x n + YH=\ a i x i satisfies 

dim Q (Q[X]/ip a (I)) = dim Q (Q[x n ]/(tp a (I)nQ[x n })). 

Corollary 3.2. Let F G Q[T], T a variable, be squarefree and r — x n +X)i=i a i x i 
with ai,...,a„_i G Z such that deg(-F) = dirriQ(Q[X]//), and F(r) G / but no 
proper factor of F(r) is in I , then I is a radical ideal. Let F — F± ■ ■ ■ F s be the 
factorization of F into irreducible factors over Q. Then I — C\ i=1 (I, Fi(r)) and 
(J, Fi(r)) is prime for i — 1, . . . , s. 

Proof. Using a linear change of variables we may assume that r — x n . Since no 
proper factor of F(r) is in / we obtain (F(x n )) = I f~l Q[# n ]. Since deg(F) = 
dimq(Q[X]/I) we have / = (xi — hi(x n ), . . . ,x n -i — h n -i(x n ), F(x n )} for suitable 
hi,..., h n -i G Q[x„]. Thus, / is radical because F is squarefree. The rest is an 
immediate consequence of Proposition 13. If I). □ 

Consequently, for the computation of the primary decomposition, we firstly verify 
whether / is already radical. Therefore we choose a generic linear form r = a\X\ + 
. . . + a„_ix„_i + x n with a±, . . . , a n -i € ^: and use a test in positive characteristic, 
similarly to section [2] 

pTestRad: We randomly choose a prime number p such that dim-p (¥ p [X]/I p ) = 
d. Let ip : F p [T] — > ¥ p [X] be defined by tp(T) = r mod p (cf. Lemma \3.blf l)) and 
(F p ) := ip-Hlp). We test tf deg(F p ) = d. 

In case of a negative result of the test there is a high probability that the ideal 
is not radical (cf. Proposition I3.1f 2^ and we compute the radical using modular 
methods. The computation of the radical is usually much more time consuming 
than the pTestRad even if the ideal is already radical. The following proposition 
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(cf. [KrLo91j . |GP07j ) is the basis for computing the radical of a zero-dimensional 
ideal. 

Proposition 3.3. Let I C Q[X] be a zero-dimensional ideal and (/,•) =10 Q[Xi] 
for i = 1, . . . , n. Moreover, let gt be the squarefree part of ft. Then the following 
holds. 

(1) VJ = I + ( 9l ,...,g n ). 

(2) //deg(/„) = dim Q (Q[X]/7) thenVI=(I,g n ). 

Proof. Part (1) of the proposition is proved in |KrLo91) . For part (2) we notice 
that if it holds deg(/ n ) = dimQ(Q[A]/7) then there exist hi, . . . , h n -i £ Q[x„] such 
that {x\ — h\, . . . ,x n -i—h n -i, /„} is a Grobner basis of J w.r.t. the lexicographical 
ordering X\ > . . . > x n . Thus, we have \fl = (x\ — hi, ... , x n -\ — h n -i,g n )- Q 

With analogous considerations as in section [21 the essential idea of the algorithm 
to compute the radical of / is as follows. Choose a set P of prime numbers, compute, 
for every p £ P, monic polynomials /{ , . . . , /„ satisfying (f> ) = I p fl F p [xj] for 
i = 1, . . . ,n and finally lift these polynomials via Chinese remainder algorithm and 
Farey rational map to (fi, . . . , /„) £ Q[a;i] x . . . x Q[x„]. 

Definition 3.4. Let (/i, . . . , /„) £ Q[xi] x . . . x Q[x n ] satisfy (fi) =1(1 Q[xi] for 
i = l,...,nE 

(1) If (A (p) , . . . , & p) ) £ W p [x x ] x ... x ¥ p [x n ] satisfies (^ p) ) = I p n F p [x 4 ] for 
i = 1, . . . , n, then the prime number p is called Zuc/cj/ for I if and only if 
deg(/i) = deg(/j ) for i = 1, . . . , n. Otherwise p is called unlucky for I. 

(2) A set P of lucky primes for / is called sufficiently large for I if and only if 
ripepP — m ax{2 • |c| 2 | c coefficient occuring in fi,. . . , f n }- 

After having computed the set FP := {(/-[ , . . . , /„ ) | p £ P} we delete the 
unlucky primes in the following way. 

deleteUnluckyPrimesRad: We define an equivalence relation on (FP,P) by 
(F(p\ P ) - (f(«), q) :^=> deg(/ l (p) ) = deg(/f g) ) /or t = 1, . . . ,n. Then the equiva- 
lence class of largest cardinality is stored in (FP, P) , the others are deleted. 

With the aid of this method we are able to choose a set of lucky primes with high 
probability. A faulty decision will be compensated by the subsequent test whether 
fi e I for i = l,...,n. 

Since we cannot predict if a given set of primes P is sufficiently large for /, we 
have to proceed by trial and error as already described in section [2j 

Algorithm [2] computes the radical of /I 11 ! 

If the pTestRad is positive then, with high probability, after a generic coordi- 
nate change it holds dimq(Q[x n ] / (I n Q[x n ])) = d. In this case it is not necessary 
to compute the radical of / and we rely on the following corollary. 



u By abuse of notation we use the same terminology as in Definition [2TT] since it is always clear 
out of context which definition we are referring to. 

The corresponding procedure is implemented in Singular in the library 
asspr imeszerodim .lib. 
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Algorithm 2 zeroRadical 



Input: I = (Gi) C Q[X] a zero-dimensional ideal generated by a Grobner basis 

Gi w.r.t. some global ordering. 
Output: G C Q[X] a Grobner basis of the radical of / w.r.t. a degree-ordering. 

choose P, a list of random primes; 

FP = 0; 

loop 

for peP do 

compute monic polynomials /,- such that (/,• ) = / p f~l F p [xj] for i = 
1, ...,n; 

FF = FPU{(/« ri ,...,/i P) )}; 
(FP,P) = deleteUnluckyPrimesRad(PP, P); 

lift (PP, P) to (/i, . . . , /„) G Q[«i] x ... x Q[x n ] by applying Chinese remainder 
and Farey rational map; 
use Gi to test if fi € I for i = 1, . . . , n; 
if /j G / for all i = 1, . . . , n then 

exit loop; 
enlarge P; 
for i = 1, . . . , n do 

compute <7i, the squarefree part of /^; 

I = I + (gi,---,9n); 

compute G C Z[X], a Q[X]-Grobner basis of I w.r.t. a degree-ordering I 1 1 
return G; 

Corollary 3.5. Let i C Q[X] &e a zero-dimensional ideal and r — x n +X)"=Ti a i x i 
with ai,...,a n -i G Z. Let P G Q[P], T be a variable, such that deg(P) = 
diniQ(Q[A"]/J) and F(r) G / fewt no proper factor of F(r) is in I. Moreover, let H 
be the squarefree part of F. Then \fl = (I,H(r)). 

Proof. The proof is a consequence of Proposition 13.31 ^2') and Corollary 13.21 □ 

Consequently we need to obtain a polynomial P G Q[T] satisfying the required 
properties of Corollary 13.21 resp. Corollary 13.51 The following lemma is helpful in 
this direction. 

Lemma 3.6. Let K be a field 1. P G Jf [P], T a variable, be monic and squarefree, 
let r = s„ + Yn=i a i x i> a i> • ■ • > «n-i G -^ suc/i i/iai deg(P) = dim.K(K[X]/I) and 
F(r) G / but no proper factor of F(r) is in I. 

(1) Let ip : K[T] ->• if [X] &e denned &y <p(T) = r. TTien y _1 (I) = (P). 

(2) Let ip : K [X] — > K [X] be defined by ip{ x i) — x i f or i = 1, . . .,H — 1 and 
■>p(x n ) = 2x n - r. Then ip(I) n K[x„] = (F(x n )). 

(3) Let A : K[X]/I —¥ K [X ]/I be the map defined by the multiplication with r, 
X(g + i) = r ■ g + i. Then F is the characteristic polynomial of A. 

Proof. (1) Since ip (P) = F(r) G i we obtain P G <£> _1 (Pj. Thus we have 

(P) = <^ _1 (i") because no proper factor of F(r) is in i. 



Here we use the procedure modStd as described in section [2] 

We substitute Q by an arbitrary field K since we also need the results of Lemma 13,61 for 
finite fields. 
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(2) It holds F(x n ) — i/j(F(r)) G i>{I) by definition of if). The assumption 
implies that no proper factor of F(x n ) is in ip(I), i.e. (F(x n )) — ip(I) H 
K[x n }. 

(3) Using the map -0 of (2) we may assume r = x n . As in the proof of 
Corollary 13.21 we obtain / = (xi — hi, ■ ■ ■ , x n ~i — h n -i, F(x n )) for suit- 
able hi,...,h n -i G K[x n ] since deg(F) = Aim.K{K[X)/I) = d. Hence, 
we may choose {1, x n , • ■ • , £^ -1 } as a basis of K[X]/I = K[x n ]/ (F(x n )), 
and obtain the polynomial F to be the characteristic polynomial of the 
multiplication with x n . 

D 

Lemma 13.61 shows that the approach of Eisenbud, Hunecke, Vasconcelos (cf. 
[EHV92] ) using (1) of the lemma, the approach of Gianni, Trager, Zacharias (cf. 
GTZ88 ) using (2) of the lemma and the approach of Monico (cf. |M02j ) using 
(3) of the remark are in principle the same. The computations for (1) resp. (2) 
require Grobner bases with respect to suitable block-orderings whereas in (3) we 
do not need a special ordering for the Grobner basis but we have to compute a 
determinant. All three algorithms are implemented in Singular. 

Remark 3.7. We can also compute the polynomial F £ Q[T] using modular meth- 
ods. For this purpose we compute F^ € F P [T] monic such that (^F^^ — ker(ip p ), 
whereat (p p : F p [T] — > ¥ P [X]/I P , f P {T) = r mod I p , for several prime numbers p 
and preserve just those F^ with deg(F^ p ') = d. Afterwards we lift the results to 
F € Q[T] by applying Chinese remainder and Farey rational map. 

Remark 3.8. If K — C is the field of complex numbers we can use the polynomial 
F of Corollary 13.21 to compute the zeros of the ideal /. The zeros of F are the 
eigenvalues of the multiplication map A defined in Lemma 13.61 Let Ai, . . . , Ad be 



the (different) eigenvalues of A then I = {~\ i=1 (I, r — Aj). Moreover, (I, r — Ai) is a 
maximal ideal in C[X] representing a zero of I for i = 1, . . . , d. 

Referring to Proposition 13. 1[ Corollary 13. 2| and the above considerations, Algo- 
rithm [3] computes the associated primes of III 

Remark 3.9. The presented versions of Algorithms [5] and [3] are just pseudo-code 
whereas their implementation in Singular is optimized. E.g., the polynomials 
fi G Fp[a:i] resp. F^ G F P [T] for p £ P are not computed repeatedly, but stored 
and reused in further iteration steps. 

Remark 3.10. Algorithm [2] resp. Algorithm [3] can easily be parallelized by comput- 
ing the polynomials f\ € F p [xi] resp. F^ £ F p [T] in parallel. Experiments indi- 
cate that the difficult and time consuming part of Algorithm [3] is the test whether 
F(r) £ I and the computation of i*i(r),. .. ,F s (r). These s + 1 computations are 
independent from each other such that they can also be verified separately in par- 
allel. 

Following the idea of one of the referees we tried to avoid the computation of F(r) 
by computing a Q[X, T]-Grobner basis of (I, T — r) w.r.t. an elimination ordering 
(eliminating X) by using modular methods (cf. section [2]) and FGLM-algorithm 
(cf. [FGLM93] ). In this case we directly compute (I, T - r) Q[X]T] nQ[T] = (F) and 



The corresponding procedures arc implemented in Singular in the library 
asspr imeszerodim .lib. 
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Algorithm 3 assPrimes 



Input: 7 C Q[X] a zero-dimensional ideal. 

Output: L = {Mi, ..., M s }, M t prime and \J1 = f)f =1 M 4 . 

compute G C Z[X], a Q[X]-Grobner basis of 7 w.r.t. a degree-ordering Pi 
compute d = dimQ(Q[J s f]/7) using G; 

choose oi, . . . , a n -i 6 Z randomly, r = aiXi + . . . + a n _ia; n _i + #„; 
if not pTestRad(c2, r, G) then 
G = zeroRadical(G); 
d=dim Q (Q[X]/<G}); 
choose P, a list of random primes; 
PP = 0; 
1 = 0; 
loop 

for p S P do 

compute P^ <E F p [T] monic such that (^p( p )) = ker(y> p ), whereat ip p : 
¥ P [T] — ► F P [X]/7 P , <^(T) = r mod 7 P 
if deg(P^) = cf then 
PP = PPU{P^}; 
if #(PP) = I then 
G = zeroRadical(G); 
d = dim Q (Q[X]/(G)); 

choose ai, . . . , a„_i G Z randomly, r = ai£i + . . . + a n _ 1 x„_ 1 + £„; 
else 

lift (PP, P) to P € Q[T] by applying Chinese remainder and Farey rational 

map; 

factorize P = Pf 1 ■ • • F" s with Pi, . . . , P s irreducible; 

compute P(r) and Pi(r), . . . , F s (r); 

if P(r) G 7 then 

if no proper factor of F(r) is in 7 then 

return {(7,Px(r)), . . . , (7, P s (r))}; 
else 

choose a non-trivial factor H of P of minimal degree such that 7f (r) G 7; 

let P^ , . . . ,Fi t correspond to 77; 

return ASSPrimes((7, F i± (r))) U . . . U ASSPRiMES((7,P; t (r))); 
enlarge P; 
I = #(7 1 P); 

may consequently omit the verification. Experiments showed that this is as time 
consuming as the presented method in Algorithm [3] 

Remark 3.11. Knowing the associated primes it is easy to compute the primary 
ideals using the method of Shimoyama and Yokoyama (cf. JSY96J ): Let Mi , . . . , M s 
be the associated primes of the zero-dimensional ideal 7 and a\, . . . ,<r s a system of 
separators, i.e. a (fc Mi and a £ Mj for j ^ i, then the saturation of 7 w.r.t. <7j is 



5 Here we use the procedure modStd as described in section [2] 
All approaches mentioned in Lemma l3,6l are applicable to verify this step. 
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rrij whereat 



the primary ideal corresponding to M,. Each oi can be chosen as n« 

rrij is an element of a Grobner basis of Mj which is not in M, . The saturation can 

be computed modularly, similarly to MOdStd and in parallel. 



4. Examples, timings and conclusion 

In this section we provide examples on which we time the algorithms modStd (cf. 
section [2| resp. assPrimes (cf. section [3|) and their parallelizations as opposed to 
the usual algorithms std resp. minAssGTSO implemented in Singular. Timings 
are conducted by using the 32-bit version of Singular 3-1-2 on an AMD Opteron 
6174 with 48 CPUs, 800 MHz each, 128 GB RAM under the Gentoo Linux operating 
system. All examples are chosen from The SymbolicData Project (cf. |G10j ). 

Remark 4.1. The parallelization of our modular algorithms is attained via multiple 
processes organized by Singular library code. Consequently a future aim is to 
enable parallelization in the kernel via multiple threads. 

We choose the following examples to emphasize the superiority of modular stan- 
dard basis computation and especially its parallelization: 

Example 4.2. Characteristic: 0, ordering: dp0, Cyclic_8.xml (cf. [BF91J - ). 

Example 4.3. Characteristic: 0, ordering: dp, Paris. iliasl3.xml (cf. [Ko La99j ). 

Example 4.4. Characteristic: 0, ordering: dp, homog. Cyclic_7.xml (cf. |BF91j V 

Example 4.5. Characteristic: 0, ordering: d:£j, Steidel_l .xml (cf. |Pf07j ). 

Tablc[T]summarizes the results where modStd* (n) denotes the parallelized version 
of the algorithm applied on n cores. In all tables, the symbol "-" indicates out of 
memory failures. All timings are given in seconds. 



Exmp. 


std 


modStd 


modStd* (4) 


modStd* (9) 


modStd* (30) 


14.21 


- 


8271 


4120 


2927 


1138 


14.31 


37734 


1159 


676 


580 


380 


14.41 


3343 


3436 


886 


408 


113 


14.51 


- 


6 


3 


3 


3 



Table 1. Total running times for computing a standard basis of 
the considered examples via std, modStd and its parallelized vari- 
ant modStd* (n) for n = 4, 9, 30. 



The basic algorithm std runs out of memory for examples 14.21 and 14.51 As 
mentioned in section [21 it is possible to parallelize the computation in several parts 
of the algorithm modStd. In many cases it turns out that the final test - the 



The procedure minAssGTZ is implemented in Singular in the library primdec.lib and com- 
putes the minimal associated prime ideals of the input. 

Degree reverse lexicographical ordering: Let X a ,XP € Mon(X). X a >d p X& :<^=> 
deg(X Q ) >deg(X^) or (deg(X a ) = deg pf 3 ) an d 31 < i < n : a n = /3„, . • • , "i-i = ft-i.c* < 
ft), where deg(X a ) = a x + . . . + «„; cf. |GP07| . 

Negative degree reverse lexicographical ordering: Let X a ,X@ £ Mon(X). X a ~>ds X& :^=> 
deg(X a ) < degiXl 3 ) or (deg(X a ) = deg pf' 3 ) an d 3 1 < i < n : a n = n ,...,CH-i = $i-i,on < 
Pi), where deg(X Q ) =ai + ... + a n ; cf. [GP07J . 
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verification whether the lifted set of polynomials includes the input and is itself 
a standard basis, see also Remark 12.51 - is a time consuming part. Therefore we 
extract the timings for the computation without the verification test in Table [21 
again in seconds. 



Exmp. 


modStd w / v . 


modStd* , (4) 

w/o v. v. / 


modStd* , (9) 

W/O V. v / 


modStd* , (30) 

w/o v. V / 


14.21 


7929 


3751 


2698 


920 


14.31 


941 


614 


552 


370 


14.41 


52 


38 


31 


36 


14.51 


6 


3 


3 


3 



Table 2. Running times for modStd and modStd* (n) with n 
4, 9, 30 without verification test. 



We consider the following examples for the computation of the associated prime 
ideals of a given zero-dimensional ideal : 



Example 4.6. Characteristic: 0, ordering: dp, Becker-Niermann.xml (cf. [D GP98] ). 
Example 4.7. Characteristic: 0, ordering: dp, FourBodyProblem.xml (cf. |BM10j ). 
Example 4.8. Characteristic: 0, ordering: dp, Reimer_5.xml (cf. |BM10j ). 
Example 4.9. Characteristic: 0, ordering: lp 20 l. ZeroDim.example_12.xml (cf. 

mm- 



Example 4.10. Characteristic: 0, ordering: dp, Cassou_l.xml (cf. [BMTO]). 

Using modular methods via the algorithm assPrimes we apply all three variants 
mentioned in section [3] 

(1) approach of Eisenbud, Hunecke, Vasconcelos (cf. [EHV92 ), 

(2) approach of Gianni, Trager, Zacharias (cf. |GTZ88j h 

(3) approach of Monico (cf. [M02] V 

We summarize the results of the timings in Table [3] and |4] where assPrimes* (n) 
denotes the parallelized version of the algorithm applied on n cores. 



Exmp. 


minAssGTZ 


as 

(1) 


sPrim 

(2) 


es 
(3) 


ass! 

(1) 


^rimes 

(2) 


*(4) 
(3) 


ass! 

(1) 


^rimes 

(2) 


*(9) 

(3) 


[____ 


- 


1 


1 





1 


1 


1 


1 


1 


1 


14.71 


- 


169 


169 


188 


104 


98 


104 


95 


100 


105 


14.81 


- 


129 


131 


230 


90 


87 


114 


76 


77 


103 


14.91 


189 


4 


5 


5 


10 


8 


8 


8 


8 


8 


14.101 


589 


35 


35 


35 


24 


23 


19 


25 


24 


25 



Table 3. Total running times for computing the associated prime 
ideals of the considered examples via minAssGTZ, assPrimes and 
its parallelized variant assPrimes* (n) for n = 4, 9. 



20 Lexicographical ordering: Let X a ,X? _ Mon(X). X a > lp X 13 
,..., ai_i = ft_i, a, > ft; cf. [GP07] . 



3 1 < i < n : ai 
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The usual algorithm minAssGTZ runs out of memory for examples I4.6[ 14.71 and 
14.81 Analogously to the modular standard basis algorithm, we also list the timings 
needed for assPrimes resp. assPrimes* (n) without the final verification step - the 
check whether F(r) £ I and the computation of -Fi(v), . . . , F a (r), see also Remark 
EZU- in Table g] 



Exmp. 


assP 
(1) 


rimes w/ 
(2) 


o vcr. 

(3) 


assPr 

(1) 


imes* , 

w/o 

(2) 


(4) 

ver. ^ / 

(3) 


assPr 

(1) 


imes* , 

w/o 

(2) 


vc,(9) 

(3) 


14.61 


1 


1 





1 








1 


1 


1 


14.71 


15 


14 


34 


7 


7 


13 


5 


5 


15 


14.81 


41 


37 


139 


39 


38 


64 


30 


26 


55 


14.91 


4 


5 


5 


9 


8 


8 


8 


8 


8 


14.101 


7 


6 


7 


5 


5 


5 


5 


4 


6 



Table 4. Running times for assPrimes and assPrimes* (n) with 
n — 4,9 without final verification step. 



Conclusion 4.11. 



(1) For the computation of Grobner bases resp. standard bases of ideals / C 
Q[X] w.r.t. global resp. local orderings modStd should be used. This is 
usually faster even without parallel computing. 

(2) The probabilistic algorithm to compute standard bases works without any 
restriction to the ordering. It is much faster than the deterministic one. 
It can be used to obtain ideas in Algebraic Geometry and other fields by 
computing several examples, similarly to computations in positive charac- 
teristic 20 years ago when computations of standard bases in characteristic 
zero have been impossible resp. too slow. 

(3) A kernel-implementation of modStd could speed up the modular part 
using the trace-algorithm of Traverso (cf. |T89j ). 

(4) An increasing number of cores used during the parallel computation of 
standard bases resp. associated primes speeds up the computation if the 
corresponding problem in positive characteristic takes some time to be com- 
puted. If the computations in positive characteristic are fast then an in- 
creasing number of cores may slow down the computations because of too 
much overhead. 

(5) In the current implementation Chinese remainder and Farey fractions are 
not parallelized. Experiments (e.g. the computation of the Grobner basis of 
Cyclic_9) show that the computations in positive characteristic need differ- 
ent time on different cores. Therefore one should apply Chinese remainder 
and Farey fractions already to partial results. 

(6) For zero-dimensional primary decomposition the modular approach is very 
efficient. This should be extended to higher-dimensional ideals. 
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